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The low-energy effective Hamiltonian of three coupled spin chains with periodic boundary con- 
ditions (spin tube) is expressed, in the limit of strong interchain coupling, in terms of XXZ chains 
coupled by biquadratic exchange interaction. A similar effective model was also proposed to describe 
the coupling of spins to orbital degrees of freedom in materials such as NaV2 0s. We investigate 
the effective model by means of bosonization and renormalization group techniques, and find that 
the generic phase diagram comprises a gapless region and gapped regions consisting of a spin liquid 
phase and various antiferromagnetic phases. We discuss the properties of the spin liquid phase, 
in particular the nature of the ground state and of the elementary excitations above it. We then 
study the effect of a magnetic field, and conclude that a strong enough magnetic field can suppress 
the dimerized phase leading to a two component Luttinger liquid. The critical exponents at the 
transition gapful-gapless are calculated and shown to be non-universal in the spin tube case or the 
generic spin orbital problem. 
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I. INTRODUCTION 

Coupled spin 1/2 chains have attracted much attention lately due to the large number of their experimental 
realizations, as well as to the variety of theoretical techniques, both analyticalj-and numerical, available to study the 
relevant models. Spin 1/2 chains, owing to the Jordan- Wigner transformation^, show properties remarkably similar 
to those of interacting one dimensional fermions with their low energy properties described by an effective Luttinger 
liquid theoryoa. Recently it has been realizedj-that these properties are drastically modified when the spin 1/2 chains 
are coupled together forming ladder systemso. In this case, in a way very similar to Haldane's spin-S problemtj, 
a gap is found to open for an even number of chains while the system remains gapless if the, number of chains is 
odd. This phenomenon has been thoroughly investigated both analyticallyQand numerically!], and corresponding, 
experimental systems were identified. Typical examples of two-chain spin 1/2 ladders exhibiting a gap ane SrC^Osu 
and Cu2{C^,Hi2N2)2ClM- On the other hand, an example of a gapless three-chain ladder is SV2 Cu^OM-, From the 
theoretical point of view, the difference between odd and even number of legsL3 or odd . an d-even spintil has been 
understood qualitatively through a generalization of the Lieb-Schultz-Mattis TheoremtSrEJIlj'til. The theorem states 
that if the ground state is unique the system will necessarily be gapless. This is indeed the case for an odd number 
of coupled spin 1/2 chains. 

More recent theoretical work has emphasized the role of boundary conditions in the transverse direction in the 
formation of spin gaps. The results quoted above - the opening of a gap for an even number of coupled chains, 
gaplessness if the number is odd - are valid for open boundary conditions (OBC). It turns out that when periodic 
boundary conditions (PBC) are imposed on the transverse direction a gap opens for both cases of even and odd 
number of chains although the underlying reasons and the natures of the gaps are different. In the case of an even 
number of coupled chains, the reason for the gap is the formation of spin singlets along the transverse direction, 
similarly to the case of chains with open boundary conditions. When- the number of coupled chains is odd a two 
fold degenerate dimerized ground state is obtained in the case of PBQi3 - in contrast to its uniqueness in the OBC - 
allowing for a gap in the spin excitations. The degeneracy of the ground state in the case of PBC can be understood 
as a consequence of the fact that PBC are frustrating for an odd number of legs. To date, no experimental system 
described by coupled antiferromagnetic spin 1/2 chains with periodic boundary conditions and an odd number of legs 
has been reported. 

When, further, leg-leg biquadratic interactions are included new states emerge. A system of two spin-1/2 spin_chains 
enters a spontaneously dimerized phase with a gapped spectrum exhibiting non-Haldane spin-liquid propertiestj. The 
elementary excitations are neither spinons nor magnons, but pairs of propagating triplet or singlet solitons connecting 
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two spontaneously dimerized ground states. Subsequently more non-Haldane spin-liquid models have been proposedS. 
Recently, these models of spin ladders with biquadratic exchange have been advocatedu as possible models for the 
formation of a spin gap in NaV20s and Na2Ti 2 Sb 2 0. 

In this paper we study the low energy physics of the three-leg ladder with periodic boundary conditions (see Fig. ^). 
This model is called the spin-tube model in the following. By analyzing in detail its effective low energy Hamiltonian 
(LEH), which consists of twaawn- equivalent coupled XXZ chains of spins and chiral degrees of freedom, in the presence 
of a biquadratic exchange£3t3, we will show that the dimerized ground state of this model falls in the universality 
class of the non-Haldane spin liquids. In the case of XXZ chains with a biquadratic coupling and for J z / J > 1, we 
find that the Ising antiferromagnetic order can compete with the dimer order, and we will describe the resulting phase 
diagram. 

The organization of the paper is as follows. In Sec.|| we introduce the spin tube model. We sketch the derivation 
of the strong interchain coupling effective Hamiltonian and discuss its relation with two spin 1/2 chains coupled by a 
biquadratic interaction. We then discuss the symmetries and recall the results already known on the two chains with 
biquadratic interactions, in particular the prediction of a spin gap. 

In Sec. QI we discuss the bosonization treatment of two non-equivalent XXZ chains coupled by a biquadratic 
interaction. We show that this Hamiltonian contains a term sin 2(p% sin 2(/>2 that is responsible for the formation of a 
spin-gap and singlet order, as well as terms cos4</>i.2 that cause antiferromagnetic order. To analyze the competitions 
of these terms, we derive renormalization group equations. Using these equations, we estimate in which part of the 
phase diagram we should expect a competition of antiferromagnetism and dimer order. We also estimate the equation 
of the phase boundary between the singlet and the Antiferromagnetic order. 

In Sec.tV we analyze in details the dimerized phase. The dimerized ground state is two-fold degenerate and is 
formed of alternating singlets of spins and pseudospins. The elementary excitations above the ground state carry a 
spin ±1/2 as well as a pseudospin ±1/2. They can be seen as formed by the introduction of a spin and a pseudospin 
in the pattern of alternating singlets. Similar excitations were obtained in Ref. |l6| and it was shown that they lead 
to response functions very different from those of a Haldane spin liquid. Numerical evidence for such excitations was 



also obtained in Ref. |20. 

In section we analyze the effect of a magnetic field strong enough to close the gap of the dimer phase . When 
both chains are magnetic, the situation is similar to the one obtained with two coupled spin laddersEircj. In the 
more physical situation, when only one of the two chains is magnetic, the ground state properties are those of a 
two component Luttinger liquid. In the case of spin-orbital models, spin and orbital degrees of freedom decouple 
completely but the presence of orbital modes should affect the specific heat. In the case of the spin-tube, we show 
that in contrast to spin-orbital models, the spin-correlation functions are affected by the presence of auxiliary gapless 
modes. We also show that for non-equivalent chains (a case that is realized in spin-orbital models) the exponent of 
the spin correlations is non-universal at the transition. Finally, in Sec. VI we will give some concluding remarks. 



II. THE SPIN-TUBE MODEL AND COUPLED XXZ CHAINS 
A. The three chain ladder with periodic boundary conditions 



We wish to study the three-chain ladder Hamiltonian, 

N 3 N 3 

H = J Si tP Si±i t p + Jj_ Si t p§i tP +i, (2-1) 

i=l p—1 i—1 p—1 

where p (resp. i) is a chain (resp. site) index, J is the coupling along the chain and Jj_ the transverse coupling. 
We impose periodic boundary conditions along the rungs by identifying Sf A = S , | t 1 . We call the resulting model the 
spin-tube (see Fig. [l|) since it can be realized by placing 3 spin 1/2 chains forming an equilateral triangle. In order 
to investigate the low-energy physics of the spin-tube, we consider the limit of strong interchain coupling (J± ^> J). 
This is the appropriate starting point that yields a good effective description of the properties of the spin tube-in the 
whole range of J±/J. (Starting from the opposite limit J± <C J, and treating J± as a perturbation, one findstj that 
it gives rise to relevant terms in the Hamiltonian. As a result, the initial Jj_ grows until it is of order of J at which 
point the weak coupling bosonization scheme is no more valid.) 

To derive the effective low-energy Hamiltonian, let us first consider the case J = 0. The system is a collection of 
independent rungs, each described by the following Hamiltonian 

H = J ± 0x.§2 + S 2 -S 3 + S1.S3) = J ± ((5x +S 2 + S 3 ) 2 - 9/4)/2. (2.2) 
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The ground state of each rung is four-fold degenerate, composed-pf two doublets of spin 1/2 excitations, corresponding 
to the left and right chirality (-/+), with energy (— 3/4Jl)E3o, and the excited states form a spin 3/2 quadruplet 
with energy 3/4 J±, 

Turning on J > allows the rungs to exchange spins. In the limit J± 3> J, the quadruplet of S=3/2 excitations 
can be neglected, and only the degenerate low energy subspace of spin 1/2 states needs to be taken into account. In 
this subspace the Hamiltonian tpa nsfpr .ms into an effective Hamiltonian with a biquadratic coupling between the spin 
and chirality degrees of freedonOc3E3, 

H *ff = | E + 2(7f 7f +1 + t?t? +1 )}, (2.3) 

i 

where S is the total spin , and n are operators exchanging left and right chiralities. The original spin operators can 
be expressed in terms of the effective spin S and the chirality r in the following way: 

Q+ - _!<?+ _|_ ln 2 PQ+T- + lnPQ+ T + 
D i,p — g J l + 3J D i T l + g J D i T i 

St P = \st + IfSfrr + 2 -fStr+ (2.4) 

The spin tube model has already been investigated numericallyll! using DMRG in the case of coupled XXX chains. 
The system was shown to exhibit a spin gap A = 0.28 J, with exponentially decaying correlation functions (— ) n {S^S^) 
and (— ) n (T§ r*). This behavior was shown to be related to the formation of a dimer order. A qualitative discussion of 
the origin of this dimer order can be found in Refj£j. In Ref. |l5| the dispersion relation of excitations having r z = 0, 1 
and S — 0, 1 as a function of momentum was obtained numerically for a system of 12 sites showing only gapped 
excitations. 



In Ref. 20, a ge neral ization of the effective Hamiltonian Eq. (2.3) to coupled anisotropic spin chains has been 



derived (see Eq. (2.10). Numerical diagonalizations were performed. A gap to S z = 1 excitations was obtained 
forO < J z / J < 1.2 with results in good agreement with those of Ref. |l5| f or J z = J. It was also shown that the ground 
state was degenerate in agreement with the dimer oxder picture of Ref. |15|. The dispersion relation of excitations having 
t z = and total spin S — or S = 1 was obtainedEl It was shown to be the bottom of a two particle continuum. The 
fundamental particle, the spinon, was conjectured to have r z = ±1/2, S ==_L/2. The spinon dispersion relation was 
obtained numerically by considering a system with an odd number of sitesE3, showing that the spinons were massive 
for all momentum. 

B. Generalized spin-tube model and XXZ chains coupled by a biquadratic exchange 



The Hamiltonian ( |2.3| ) is part of the class of the Hamiltonians consisting of two non-equivalent coupled XXZ chains 
in the presence of a biquadratic exchange, 

H = H +H B , (2.5) 

where 

i a=l,2 

Hb = ^/l(SfjSij +1 + Sf^Sf^ + Ai£f )i £f >i+1 )(52 > jS2 i j +1 + ^2,1^2,1+1 + ^^2,^2^+1) ■ (2-7) 

i 

In the case of the spin-tube, Su corresponds to spin S%, and S^,* is associated w ith the chiral degrees of freedom n. 
Another way of writing a class of Hamiltonian generalizing the spin tube model ( |2.3[ ) is: 

N 

H = ^[u + j(S+Sr +1 + S rS+ +1 ) + J z SfS? +1 ] [v + a(T+ T r +1 + T r T + +1 ) + frf rf +1 ]. (2.8) 

8=1 

The effective Hamiltonian for the spin tube is obtained for, 

u = 7 = J/6 J z = J/3 

v = l a = l J' z =0. (2.9) 
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For a tube made of XXZ chains, one has, instead^ 



u = 7 = J/6 J z = J/3A 

v = l a = l J' z =0. (2.10) 



The parameters J z /l and J^/a in Eq. fl2.8|) measure the XXZ anisotropy for spin and chirality, respectively. When 
both of them are equal to 1, the Hamiltonian Eq. ( [2.8] ) is SU{2) x 5/7(2) symmetric. For u = v, a = ~f, J z = J' z the 
two chains are equivalent, and can be parameterized as: 

u = « = |, a = 7=1, J z = j z = 2q. (2.11) 



Hamiltonians of the type Eq. ( |2.8|) can be mapped onto Hamiltonians of the type Eq. ( [2.5| ) . The correspondence 
is given by: 

Ji = 2ua J2 = 2«7 

Jf = v J z J 2 2 = uJ' z 

J V 

27 2a 

A = 4a 7 (2.12) 



with the identification Si = S , S2 = t. Since writing the Hamiltonian in the form (2.5) is less restrictive than in 
the form (|2.8[) (i.e. the former includes the case A a ^ y^), we focus on the Hamiltonian (2J)). Hamiltonians of the 
type Eq. (|2.5|) are also encountered in a different context than the spin-tube model. In particular, a Hamiltonian p£ 
the type (2.5) has-been proposed by Mostovoy and Khomskii as a model for the spin gap formation in the NaV205t3 
and Na 2 Ti 2 Sb20c3 compounds . In that case, the Si spins correspond to the real spin of the system whereas the S2 
spins are pseudospins associated with orbital degrees of freedom. These spin-orbital models can be derived from a 
multiband Hubbard model. The derivation is reviewed for instance in Ref. 2^. 

let us discuss first the case J\ = J2 = Jf = Jf = J, A x = A 2 = 1. The Hamiltonian describes two coupled 
Heisenberg chains with a biquadratic coupling preserving the SU(2) symmetry. Actually, the full symmetry group is 
larger than SU(2). One has: 

[H,Si,tot.}=0 [H, S 2 ,tot.] = (2.13) 



and the full symmetry group is therefore SU (2) x SU(2) rather than the SU(2) symmetry that follows from [H, Si,, 
$j,tot.] =0. As a result the spectrum consists of SU(2) x SU(2) — SO (4) multipletaia. For J — j, the Hamiltonian 



( [2.5| ) has been shown to have an even larger St/ (4) symmetry and to reduce tjQ-iaa integrable St/ (4) spin chainE3. 
The spectrum of the SU(4) spin chain has been obtained by, the Bethe AnsatzEScEl, and the correlations functions 



have been obtained by non-abelian bosonization techniquesE-J, identifying the low energy effective theory describing 
the spin chain as the St/(4)i WZNW model. JEhe integrable St/(4) spin-chain has also been intensively studied 
numerically using Quantum-Monte Carlo (QMC)E3 or Density Matrix Renormalization Group (DMRG) and Lanczos 
Exact Diagonalization (ED)E3 in the context of the spin-orbital models. The ground state energji-and excited state 
energy were obtained in good agreemeriJtEp with analytical calculation using the Bethe Ansatzc3. The numerical 
calculation of the correlation functionsE3H reproduces results of the the continuum field theoretical treatmentEi At 
J A/4 perturbations are generated which lower the St/(4) symmetry to St/(2) x SU(2) and render the Hamiltonian 
non-integrable. These perturbations have been recently studiedLj by describing the chain away from the integrable 
point as a perturbed St7(4)i WZW model. It was found that for J < A/4, a gapless phase is obtained while for 
J > j, a gap is formed. A different field theoretical treatment in the limit A <C J also leads to the appearancejaf 
a gapped phase. For J = 3A/4, the ground state wave-function could be obtained exactly in matrix product formliU, 
of singlet states along the legs of the ladder. This picture is in good agreement with the predictions of the field 
theoretical treatment. The dependence of the gap on the coupling for the range < A/J < 4 was obtained by DMRG 
calculationsc3i2a. It was found that for A/J -C 1, the gap increases proportionally to A in agreement with the peak 
coupling bosonization treatment and vanishes for A/J = 1/4, which is the SU(4) symmetric point as predictedc3E3. 
The DMRG calculations of Ref. Eq showed however a power law gap opening. Such power law gap opening can only 
be explained by the presence of a relevant operator in the continuum description. However, no such operator was 
obtained in the bosonization treatmentE-i Moreover, if a relevant operator was present in the continuum theory, the 
absence of a gap at A/J = 4 would be the result of the coefficient of this operator vanishing precisely at A/ J = 4. But 
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then, in contrast to what is observed in numerical calculations, a gap would also obtain for A/ J > 1/4. A solution 
to this puzzle has been suggested recently^! In Ref. [2(| the gap has been calculated by assuming that the first 
excited state was in the subspace (Si, S2) = (1,1)- This assumption was shown to ha,incorrect in the gapped phase 
in which the first excited state lies in the subspace (Si, S2) = (1, 0). When corrected!^!, a slow increase with 4 — A/ J 
is found, compatible with an exponential gap -opening. The correlation functional do not show incommensuration, 
in agreement with the field theoretic approacto. 

Taking J\ — Jf ^ J2 ffn^f an d Ai = A2 = 1 in ( fE5| ) preserves the SU(2) x SU(2) symmetry. This case has been 
investigated numerically^ with J\ = 0.7 J2 — 0.3A/4. It was shown that a gapless-gapped phase transition obtained 
as J2/A was increased This work was followed by analytical investigation based on the perturbed SU(A)i WZW 
continuum theoryEjO. The analytical investigations established the existence for a given A of an extended gapless 
region in the plane Ji J 2 that contains the line J\ — J2 < A/4 previously discussed. 

In the present work, we consider the general case of J* ^ J a and A Q 7^ 1. This case includes in particular the spin 
tube model. We will focus on the regime A <C Ji,2 and apply methods similar to those of Ref. pl| 



III. PHASE DIAGRAM 



In this section, we derive the phase diagram of two XXZ s pin chains weakly coupled by a biquadratic exchange. 
We first recall the bosonization of a spin chain in Sec. Ill A, This section can be skipped by readers familiar with 

This allows us to 
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bosonization. Then, we discuss the bosonization of the two coupled chains system in Sec 

derive Renormalization Group (RG) equations for the coupling constants of the problem. Finally, in Sec. fIID|, we 
discuss the phase diagram deduced from the analysis of RG equations. 



A. bosonization of a single XXZ chain 

In this section, we recall briefly the derivation of the bosonized Hamiltonian and smp_aaexators that describe an 
isolated XXZ chain. We follow the well known abelian bosonization procedure for spinsEfEjO'LJ. The XXZ spin chain 
is described by the Hamiltonian: 

ffxxz = j£>fSf +1 + SfSf +1 ) + J 2 £ SiS? + i (3-1) 

i i 

The XXZ spin chain Hamiltonian is first transformed into an interacting fermionic system on the lattice by expressing 
the spin operators S + , S~ , S z in terms of fermion operators (V, a, using the Jordan- Wigner transformation^: 

i-l 

S+ = (-)*<»} cob [7r5>t 0j ] (3.2) 

3=0 

Si - (-Y COS f 7T J2 a ] a A Q i ( 3 - 3 ) 

St=a\a t - 1 -. (3.4) 

This transformation turns the XXZ Hamiltonian into a model of spinless fermions with nearest neighbor interaction 
described by the Hamiltonian, 

Hxxz = -^r^2{a\a l+1 + a\ +1 Oi) + J z ^(ajctj - x)(af+i a i+i - h. (3.5) 



For J z = 0, the Hamiltonian (3.5) describes non-interacting fermions and is easily diagonalized. To proceed, we 
restrict our attention to the low-energy sector of the theory, captured by the continuum theory. Introduce the left 
(right) chiral fermion fields iPl(%) (V'-R( a; )) containing momenta close to the Fermi points kp = i^S (x = na with a 
the lattice spacing), 

^= ^e'T^H + e-'T^fna). (3.6) 
Ja 
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The continuum Fermi fields are then reexpressed in terms of bosonic as follows: 

e %(9-<t>){x) 



\l2ixa 

e i(e+4>)(x) 



l2ixa 



where the pair of conjugate fields, II, <f> satisfy the following commutation relation: 

[$(x),II(a;')] = i5{x-x'), 

and the field 6, dual to c6, is defined as: 



The spin operators Eqs. ( p.Sj )— ( 3.4 ) , can be expressed in the continuum limit as: 



S + (x) = 



S,„ 



i0(x) 

>2ira 



cos 2c 



d. 



■ cos2c6, 



(3.7) 



(3.8) 



(3.9) 



(3.10) 



Introducing normal ordering with respect to the fermion vacuum, one has: 

ox qx , qV qV _i_ A Q z Q z / Q x Q x j_ oV cV i A Q z Q z \ _l 

, . ca; Q£ I of CV i A Q z Q z 



(3.11) 



where : . . . : indicates normal ordering. The average_/l = (SfSf +1 + Sf Sf +1 + ASfS* +1 ) is the ground state energy 
of the chain and can be found by the Bethe AnsatzEl. However, as long as one is only interested in the correlation 
functions of the single chain, one can simply drop this contribution and focus on the normal ordered terms that 
describe the excitations abov e the ground state. Our task is thus to derive a bosoaiged expression of the normal 
ordered product in Eq. (3.11). Some care is needed in order to obtain correct resultscireHl, and one finds: 



. Q z Q z i Q z Q z 
■ D i D i+1 ~r D i 



- [(nil) 2 + (O x( f>f 

IX 

\{d x <P) 2 + 

ir z (TTa) z 



TV 

' sin 2(j> 



sin 20 

2 cos z 



(2™) 



irrelevant terms 



(3.12) 



The oscillat ing terms in Eq. (3.12) are dropped from the Hamiltonian after integration over x. The Hamiltonian 
Hxxz, Eq. ( |3.5| ) then becomes: 



H 



xxz 



dx 
2^ 



uK(irK) 2 



K 



where for J z <C J, 



aJ 1 



(d x <t>? 



4J 2 



28 



(2™) 



da;cos(40), 



ttJ 



1/2 



K = ^1 H 
S = J 7 a. 



irJ 



-1/2 



(3.13) 



(3.14) 

(3.15) 
(3.16) 



Thus, the bosonized form-Qf-J/xxz reduces to a sine-Gordon Hamiltonian, .where the cosine terms come from the 
intrachain Umklapp processc^la. The renormalization group treatment showsc3 that in the vicinity of the XY point, 
the cosine terms are irrelevant, so that asymptotic properties are described by a free scalar field with renormalized 
u* , K* . Since the XXZ chain is integrable, it can be shown that the gapless spectrum extends to the whole region 
\J Z \ < J— Moreover, it is possible to obtain an analytic expression for the renormalized u* , K* from the exact 
solutionM. One finds: 



G 



(3.17, 

2 arccos ■ 



^ = n 2„.* ( 3 - 18 ) 



2—4 arccos ^j- 



The isotropic point J z — J (Antiferromagnetic Heisenberg model) corresponds, in the bosonization description, to 
K* = 1/2 and 5* = 0. At this point cos(4<£4 is marginally irrelevant. For J z > J, the exact solution shows that a 
gap opens in the excitations of the systemucj, and that an Ising order of the spins along the z axis is obtained. This 
result can also be found from a bosonization procedure valid in the vicinity of the isotropic point. 



B. Two coupled XXZ chains with a biquadratic exchange 



We now proceed to derive, using the results reviewed in the previous subs ecti on, the bosonized Hamiltonian of two 
non-equivalent coupled XXZ chains. To bosonize the Hamiltonian Hq, Eq. (2.6) describing two decou pled XXZ chain 
we introduce two pairs of dual fields (one pair for each chain) 8 a , <p a (a = 1,2) as defined in Sec. Ill A. The bosonized 
form of the Hamiltonian Hq is then : 



a=l,2 w 



dx 
2^ 



u a K a {irTi. a ) 



2Sg 

(27ra) 2 



dx cos 40 Q 



(3.19) 



with the fields <pi and 2 having a prori different velocities and Luttinger couplings, 



K* 



2 arccos 
1 



J" 



2 — — arccos -y 1 

7T J a 



S a = J z a a. 



(3.20) 



The bosonization formulas for the spins (3.10) are unchanged except for the obvious introduction of a chain index. 

In order to have the full bosonized Hamiltonian, we now have to derive the bosonized form of the biquadratic 
exchange (2.7). The first step is to normal order using Eq. (3.11). This step is important since it leads to non trivial 
contributions to the quadratic part of the Hamiltonian. We introduce the terms: 



A, 



qV qy I a qz qz 



+ l,al 



(3.21) 



where a = 1, 2 that measure the strength of these corrections. Having expressed the Hamiltonian in a normal ordered 
form, we can apply Eq. ( 3.12j ) to obtain the bosonized expression. The oscillating terms of Eq. (3.12) give rise to an 
interchain coupling term, the strength of which is given by: 



g = 2Xa (1 + 2Ai/tt) (1 + 2A 2 /tt) . 
Collecting everything, one obtains the following bosonized biquadratic exchange: 

2g /".._,._, f dx 



(3.22) 



Hp 



(27ra) 2 



dx sin 20i sin 2<p2 + A\ 



2- 



[(7rn 2 ) 2 + (l + 2A 2 /^)(a x 2 ) 



A, 



dx 



(l + 2Ai/7r)(a^i) 5 



irrelevant terms 



(3.23) 



The terms in A a form a mean-field like interchain interaction . For Jx,2 ~/~ 0> these terms merely produce a renormal- 
ization of the velocities and Luttinger liquid exponents with respect to the decoupled chains. However, in the case 
where J\ = or J 2 = 0, these terms become crucial since they give a finite velocity to the <pi (resp. </> 2 excitations) 
excitations. Such a case is realized in the spin-tube problem where the exchange constant of the pseudospins r is 
zero. We see that the mean-field like contribution of the spin fluctuations contributes in that system to provide an 
exchange constant and thus a finite velocity to the pseudospin excitations. Of course, in such case, the bosonization 
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procedure is not really justified since interactions are of the order of magnitude of the bandwidth of the spin or 
pseudospin excitations. However, it is usual in quasi-one dimensional systems to have a continuity between the weak 
and the strong coupling regime. Moreover, a mean field theoretical treatment in the XY limit (J-f = Jf = Af = Af) 
leads to similar results to the bosonization treatment. Details can be found in App. |a|. We will therefore assume 
that although not fully justified in the spin tube case, bosonization nevertheless leads to qualitatively correct results 
concerning the phases of the system and the overall behavior of correlations in these various phases . A quantitative 
treatment (in particular of the phase boundaries) requires numerical simulations that are beyond the scope of this 
paper. 

The present treatment shows us that in weak coupling the two chains are described by a bosonized Hamiltonian: 



2tt 

a = l,2 



u a K a (TTlI a ) 2 + ^-{d x (j) a ) 2 - """" cos40 a 
K a {2Tra) z 



2.9 



(2na 



c?x sin 20i sin 202 (3.24) 



We shall analyze the model perturbatively, with results valid up to a given value of A: it is known that in the 
strong coupling regime of the two chains with biquadratic exchange there is a special value of the interchain coupling 
at which the model has SU(4) symmetry. At this special point, the model develops a gapless phase described by a 
SU(A)i WZNW model. Such an effect is non-perturbative: the resulting critical point has a confermal anomaly c = 3 
whereas the original unperturbed model has c = 2. This implies by Zamolodchikov's c-theoremc£l that the transition 
to the SU(4)i cannot be predicted by a RG calculation that always leads to a decrease of c. Beyond this special value 
of A , the weak coupling theory would lead to incorrect predictions and an alternative approach such as the one of 
Lecheminant and Azaria is needed. On the other hand, it is expected to give qualitatively correct predictions when 
the coupling is smaller than the critical value. In the remainder of the paper we shall thus work in the weak coupling 
regime where the weak coupling theory is valid. In the following section, we will discuss the RG treatment of the 
weak coupling model. 



C. Renormalization group equations 



In the following analysis, we will neglect the velocity difference between chains 1 and 2 since usually velocity 
differences do not play an important role in the derivation of the phase diagram by RG techniques. However, 
for the sake of completeness, we have-given in App. [§| the RG equations for non-equal velocities derived using a 
momentum shell integration tcchniqud£D. When velocity d iffere nc es are neglected, the renormalization group equations 
for K a ,8 a ,g can be easily derived from the Hamiltonian ( 3.1S )-( 3.23 ) using Operator Product Expansions (OPE)cfl. 



The renormalization group equations for 5 a , g neglecting velocity differences are: 

d ( 8i 



TTU 
2 

TTU 



4Ki) 


8 L _ 


9 2 




TTU 


8tt 2 u 2 


AK 2 ) 


82__ 


9 2 


TTU 


8tt 2 u 2 



dl 

d_ ( 8 
Jl 

dl \TTUJ TTU 2TT Z U 



while the renormalization group equations for K\ , K 2 are: 



±(±) = (h.y + i(± 

dl \Ki J y ttu J 8 V ttu / 



We now proceed to deduce the weak coupling phase diagram of the model. Eq. ( 3.25| ) indicates that g is a relevant 
variable when (Ki + K2) < 2, while the variables 61,62 become relevant when respectively, Ki < 1/2, K2 < 1/2. 
There are therefore four cases to distinguish. In the first case: Ki + K2 > 2, there are no relevant operators and the 
system is in a gapless state. In the second case: Ki + K 2 < 2, Ki, K 2 > 1/2, there is a single relevant operator; for 
Ki > 1/2, K 2 < 1/2 (or equivalently Ki < 1/2, K 2 > 1/2), there are two, while for Ki,K 2 < 1/2 there are three 
relevant operators. These four different regions are shown in Fig. @. 
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In the presence of relevant operators, RG equations cease to be valid as soon as the dimensionless coupling become 
O(l), and we have to determine the nature of the strong coupling fixed points in order to predict the phase diagram. 
We see that K\ , K2 are driven to zero by the flow of the RG when there is a relevant operator so that the fields 
become classical. As a result, fa, fa are lockfitLat average values (fa), (fa) that minimize the ground state energy and 
a gap opens in the excitations of these fieldsca. When these fields are locked, it is also known that the exponentials 
of the dual field have exponential decayed. 

Let us consider first the case with g the only relevant operator. Then, the minimization of the ground state energy 
requires svn2(fa) sin 2(02) = —I, i.e. (fa) = —{fa) — ±j ■ It is then clear that cos20i j 2 as well as e 16 , which are 
the staggered component of S z and S + , will display an exponential decay. It can also be shown that the uniform 
component of the spins also present an exponential decay. As a result we have a spin-liquid phase that presents 
only short range order in all its spin correlations. As we shall see, this spin-liquid phase is a dimcrized phase whose 
properties are discussed extensively in Sec. IV. In the case 3 = (decoupled chains) the analysis is even simpler. Then 
the two chains remain decoupled, and we are left with the analysis of the usual sine-Gordon model. It is then known 
that when 8 a is relevant, the field 4> a is locked. The analysis of the resulting strong coupling fixed point can be found 
for instance in Ref. ^|. It is found that the strong coupling fixed point is associated with the Ising Antiferromagnetic 
phase of the XXZ chain at J z > J in which the staggered component of S z has a non zero expectation value in the 
ground state. 

W hen q and at least one of the 8 a are relevant, there are two possible candidates for the ground state. Considering 
Eq. ( |3.25 ), we see that the effect of the biquadratic interchain interaction is to reduce the effective <5i,2- Physically, this 
means that the tendency to form singlets competes with the tendency to form an Ising antiferromagnet. Two scenarios 
are possible. One is that there is a well defined phase boundary between a pure spin-liquid state and a pure Ising 
antiferromagnet state. The second scenario is that there is a crossover between the two pure states as A/ J-f 2 is varied. 
In such case, increasing A would lead to a gradual disappearance of antiferromagnetic order leaving a purely singlet 
state as A — > oo. Since in both phases it is the same field that orders, there is a priori no reason to exclude a mixed 
spin-liquid antiferromagnet order. Thus, the first scenario appears extremely unlikely. A numerical investigation of 
the crossover could be very interesting as a toy model of a crossover fro m spin -liquid to antiferromagnetism. It is 
interesting to remark that if K\ = K2 = 0, and S\ = 82 = 8, the equations (3.25) can be integrated analytically. Two 
phases are obtained, separated by a line g = 2^/28. In the first one, g — > +00 and 8 — > —00 which corresponds to 
a ground state with singlet order. In the second one, g — > and 8 — > +00, which corresponds to antiferromagnetic 
order. The RG equations cease to be valid for Si/(iru) ~ 1 or g/(iru) ~ 1. If when this scale is reached g and 8 are 
of the same order of magnitude, there is a possibility of obtaining a mixing of antiferromagnetism and dimer order. 
Note that even on the line g = 2\[28, there is a finite correlation length. This is a further evidence for a progressive 
crossover from dominant Antiferromagnetic order to dominant dimer order. 

It is also possible to give a purely classical treatment for K\ = K2 = by simply minimizing the ground state 
energy with respect to (fa) and (fa). In the case S 2 > Si, one finds that there are three different regimes. For g > 4<5 2 
one obtains (fa) = —(fa) = ±f corresponding to a purely dimerized phase. For 4^/8162 < g < 4^2, one obtains 
(fa) — ±j and sin 2 (fa) = =Fg/(4<52). This corresponds to persistence of dimerization in the chain with the smallest 
tendency to antiferromagnetic order, whereas the chain with the strongest tendency to antiferromagnetism is found in 
state with mixed dimer and antiferromagnetic order. The antiferromagnetic order parameter in that chain, cos 2 (fa), 
then assumes the value ±yl — (g/(4<52)) 2 . Finally in the region g < 4y/SiS2, both chains display antiferromagnetism 
with (fa,2) = ±7r/2. These results are summarized on figure^. For 81 = 62, the region with mixed antiferromagnetic 
and dimer order shrinks to a single point. It would be interesting to see how quantum fluctuations affect the present 
picture and in particular determine whether sharp transitions are preserved or if they evolve into crossovers. 



D. Phase diagram in zero external magnetic field 



In this section, we try to estimate the position of the crossover between the spin-liquid and the antiferromagnet. 
This can be done roughly by comparing the correlation lengths in the antiferromagnet and in the spin-liquid phase. 
Using the RG equations, for small g, S±, 82, we can neglect the renormalization of Ki, K 2 - This leads to: 

g(l)=g(0)e^- K ^ 
S x (l) = 5i(0)eV- AKl » 

8 2 (l) = 8 2 (0)e^ iK ^ (3.27) 

when any of these quantities become of the order of the energy cutoff ttvf/ci, the RG equations cease to be valid and 
the phase that is obtained is determined by minimizing the ground state energy. 
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For g, the strong coupling is obtained at a length scale: 



V a 9 J 



2-Ki -K, 



This is the correlation length of spin fluctuations in the spin-liquid phase, 
respectively at length scale: 



Laf,i = a 
Laf,2 = a 



ttvf 

■KV F \ 



(3.28) 

For 61,62 the strong coupling is obtained 

(3.29) 
(3.30) 



ft is clear that the shortest length corresponds to the first operator to attain strong coupling. Therefore, the phase 
that is obtained is the one with the shortest correlation l ength . For K\ = K2 = 0, this is in agreement up to a 
constant with the criterion derived from the RG equations (3.25)-(3.26). The comparison of correlation length allows 
to draw a rough phase boundary between the antiferromagnet and the dimerized phase that could also be obtained by 
numerically integrating the RG equations (3.25)-(3.26) starting from weak coupling and any K\,K<i. The equation 
of the phase boundary is in the case of equivalent chains 6\ = 6% = 6, K\ = K2 = K : 



fJ 



ttvf ( \6\a 
a \ttvf 



(3.31) 



Let us note that in the isotropic spin-tube case, the operators causing antiferromagnetic order are (marginally) 
irrelevant so that there is only singlet order. However, in the case of an anisotropic spin tube (3 coupled XXZ spin 
chain with J z < J), such competition becomes possible. Completely decoupled chains exhibit for J z < J an Ising 
antiferromagnetic phase. Introducing a strong enough biquadratic interchain coupling favors on the other hand a 
spin liquid phase. The competition of the two should produce a crossover from the Ising Antiferromagnet to the spin 
liquid of the type discussed in the preceding section. 



IV. SPIN LIQUID PHASE 



In this section, we discuss the properties of the spin liquid phase. In the case of equivalent chains, further progress 
can be made by using symmetric and antisymmetric modes allowing in particular a refermionization of the problem 
and the calculation of some correlation functionsE3. This will be discussed in Sec. IV A. In the general case with 



inequivalent chains, such decoupling is no longer possible. However, it is still possible to present a sim ple semiclassical 
picture of the nature of excitations above the ground state. This will be the subject of section IV B . 

We will focus on region | < K\ = K2 < 2, K\ + K2 < 2, in which the only relevant operator is the biquadratic 
exchange, sin20i sin 2</>2- The ground state shows long-range order of the fields 4>i and (f>2- The expectation values of 
the ordered fields are: 



TT 

= ± 4' 



7T 

= T 4' 



(4.1) 



and as a result, (sin2</>i.2 
and pseudospin variables, i.e 



^ 0. Using Eq 

( 1)' (S a iSai 



(3.12), this implies that a dimerized order develops both in spin variables 
(a = 1, 2). In parallel with that, we have (cos2(/>i i 2) = 0, so that 
by Eq. (3.1C ) (— = 0, and the correlation functions (— (S^S^A — + as \i — j\ — * 00. It is also well known 



\i-j\-> 00. 

that when the fields <p a are ordered, the correlation functions of the disor der o perators e lpe ° decay exponentially at 
large distances. Using again the bosonization formulas for the spins, Eq. ( 3.10| ), this implies an exponential decay of 
all correlation functions : (— p 1- ^ (S£ S® ) and (SfSj) where a — x,y,z. Therefore, the dimerized phase appears as a 
spin liquid state formed of singlets of spins-on both chain 1 and chain 2. Such a conclusion was reached previously p. 
a numerical investigation of the spin tubecJ and by considering the equivalent isotropic chains at the solvable pointtll 
A = 3 J/4 where the ground state wavefunction can be obtained exactly in Matrix Product Form. Our results show 
that such mechanism of spin liquid ground state formation does not require SU(2) symmetry. 

This mechanism is somewhat reminiscent of the spin-Peierls trapsitiono, the pseudospins playing here the role of 
the phonons. This is the essence of the Mostovoy-Khomskii modellia for the "spin-Peierls" transition at T c = 34K in 
NaV 2 5 . 
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In the case of the spin-tube the same picture of the ground state obtains, with the ground state of the spin-tube 
formed by singlet of spins on even bonds and singlet of chiralities on odd bonds or by singlet of spins on odd bonds 
and singlet of chiralities on even bonds (see Fig. 0) . 

For the moment, we have only been able to discuss the nature of the ground state. However, it is also important 
to discuss the nature of the excitations as well as the various correlation functions. In order to do that, it is worth to 
restrict first to the simple case of two equivalent chains, in which the physical picture is the clearest. 



A. Equivalent chains 

When the two chains are equivalent, we introduce the fieldsEl: 



»1-|"W/V2, <t>a = (4>1 - 02)/V2 



and their conjugate fields: 

n s = + n 2 )/V2, n a = (n a - u 2 )/V2, 

so that the total Hamiltonian can be completely decoupled into the symmetric and antisymmetric parts, 

H = H s + H a , 



(4.2) 
(4.3) 
(4.4) 



dx r 

27 



K 



9 



(27ra) 2 



dx cos 



(4.5) 



dx 
2^ 



uK(TTU a ) 2 + ~{d x <i>af 
A 



(2na) 2 



dx cos 



(4.6) 



where the magnetic field couples only to <f> s , and only the most relevant operators have been taken into account. The 
elementary excitations can be discussed in terms of solitons of two decoupled sine-Gordon models. The solitons of the 
Hamiltonian H s carry a total spin m = ±1 whereas those of Hamiltonian H a carry a spin as they do not couple to 
the magnetic field. These solitons are represented in Fig. M. It is also convenient to use the canonical transformation 
cf> a .s = s , a /\/2, n s . a = y/2Il S:a followed by a refermionization. Introducing the fermion operators: 



^L.u{ x ) 

where v = a, s, one can finally rewrite H a s in the form: 



e l{6v-<i>v) 

\j2ita 



'2na 



iv I dx(tlr R v d x ipR tU — ip T L v d x ipL,v) — fJ-v / dx{ty\ v ty L , v + T Li $r, v ) + A / dxp v (x) 



where p v (x) — rp L v ^l,v + Wr v ^R,v and v — a,s. The couplings are given by: 

v = 2uK 

A = iru ( K 

\4K 



(4.7) 



(4.8) 



47ra 



= Ma = -Ps 



(4.9) 



At the isotropic point (K = 1/2), one has A = 0, so that H u is a free fermion HamiltonianEl. Similarly to the 
spin ladder, where [Stot.,H] = 0, the excitation spectrum can be split into a singlet and a triplet with spin —1,0,1. 
However, in contrast to the spin ladder case, the triplet and the singlet here have the same mass. This is the signature 
of a larger symmetry group, SU{2) x SU(2) ~ 50(4). The-jjojpelation functions can be obtained from mapping the 
free fermion Hamiltonian onto two non-critical Ising modelEirBlla exhibiting the SU{2) x Z2 symmetry. Remarkably, 
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although thepSystem has a spin gap, these correlation functions are very different from those of a spin-1 chained or 
a spin laddeil£3. In the chain with biquadratic exchange, the response functions do not show any particle-like delta 
function peak in their imaginary paj-t, but only a two particle continuunJla even in the vicinity of q = n. This is to 
be contrasted with the spin ladderEj which shows a delta function peak associated with a single particle excitation at 
q = 7T. The two chains with biquadratic interactions thus form a "non-Haldane" spin liquid. 

Away from the isotropic point, the Hamiltonian can still be refermionized but the fermions (solitons) have in- 
teractions which preclude a mapping on a non-critical Ising model. However, in the anisotropic case, when the 
antiferromagnetic intrachain interaction is irrelevant one has K > 1/2. This implies that no bound states of solitons 
can be formed since there is no coherent propagation of two solitonaij, hence the absence of a single particle peak. A 
calculation of correlation functions is in principle possible using as form factors approach, but this is far beyond the 
scope of the present paper. Since H a and H s remain decoupled, the excitations having m — ±1 do not interact with 
the excitations having m = 0. This can be understood as a consequence of the U(l) x U(l) symmetry of the problem 
and the resulting separate conservation of Sf and S 2 . 



B. Non equivalent chains 



In the case of two non-equivalent chains, one cannot decouple the Hamiltonian into two sine-Gordon Hamiltonians. 
Therefore, determining the nature and the quantum numbers of elementary excitations is not as straightforward. The 
preceding canonical transformation leads to the following Hamiltonian: 



H 



dx 
2^ 



(u\Ki + u 2 K 2 ) 



u x K x - u 2 K 2 2 

H 7T ll s ll a 



[(ttU s ) 2 + (7rn a ) 2 ] 
ut/Ki - u 2 /K 2 



(u 1 /K 1 +u 2 /K 2 ) 



d x 4> s d x (t> a 



■Mi 



(27ra) : 



[(O^s) 2 + (d x <f> a ) 2 } 

J g?x[cos V8(f> s — COS V8(f> a } 



Shifting (f> a — ► 4> a + -^=, we see that the resulting Hamiltonian has a Z 2 invariance under (n a ,</> a ) 
rescaling (f> a — ij> a /^/2, tfi s — s /y2, and a refermionization brings the Hamiltonian to the form: 



(n s 



(4.10) 
6 S ). A 



H 



dx 



r—s,a 



+ 9{PR.s{ X )PrA X ) + PL,s(x)pL,a(x)) +91 ^ PR,r(x)pL,r'(x) + <?2 ^ P 'rA X ) 'P 'L ,r( X ) 



(4.11) 



As a result, there is now an interaction between the excitations of spin S z — (the a fermions) and the excitations 
of spin S z — ±1 (the s fermions). We have: 



v = — 



ui K x 



1 



g = 7T 
91 



7T 



ui ( Ki 
1 



.92 



m 



7T 



Ui 



4Ki 
1 



4Xi 
1 

-Ki 



- u 2 yK 2 + 
u 2 (k 2 + 
u 2 ^K 2 — 
u 2 \K 2 - 



1 



4JsT 2 
1 

4^ 
1 

AK~ 2 
1 

4^2" 



9 
Ana 



(4.12) 

(4.13) 

(4.14) 

(4.15) 
(4.16) 



The fermionic version of the model is a generalization of the massive Thirring model. The fields carry spin and the 
interactions break spin rotation symmetry. It can be checked that the interaction of a with s fermions disappears only 
for Ui — u 2 , Ki = K 2 . Not much can be said of the elementary excitations of Hamiltonian (4.11) due to the absence 
of an exact solution. In order to gain some insight into the elementary excitations of the dimerized state in the case of 
no n-equ ivalent chains, we resort to semi-classical approximations. Clearly, we can search for a semi-classical minimum 
of ( |4.1CD with either (f> s = or <j> a — Tr/y/8. This corresponds to a single soliton in the system. In the general case, 
such excitations are associated with a single fermion, either a or s. In order to obtain a physical picture for this type 
of elementary excitations, let us calculate the average "magnetization" for the fields (j>\ and 4> 2 when there is a soliton 
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connecting the two degenerate dimerized ground-states. In the case 4> s = 0, 4>i decreases from tt/A to — tt/A while 4>2 
increases from — tt/A to tt/A. We immediately get: 



1 f°° 1 1 

mi = — / d x 4> x = — ^i(+oo)-^i(-oo) =-- (4.17) 

7T . _^ TT 2 



771 2 



i r°° i i 

— / 3 X 2 = ^2(+00)-^2(-00) = -TO! = -. (4.18) 

TT /.^ TT 2 



Another solution is obtained by reversing the sign of (j) a leading to mi = —7772 = 1 /2. The case <j) a — tt/^/8 corresponds 
to mi = fni = ±1/2. In each case, the elementary excitations are formed by breaking a singlet on neighboring sites 
on each chain. Such objects then propagate coherently. In the case of the spin tube, this corresponds to having one 
unpaired spin associated with one unpaired chirality pseudospin. Such an excitation has r 2 = ±1/2, S z = ±1/2. It is 
the spinon of Ref. ^0|. Thus the elementary excitations of the model can be easily visualized as an unpaired spin and 
an unpaired chirality forming a triplet or a singlet diagonal bond (see Fig. |7|). This is the soliton in the dimer order. 
This picture generalizes from isotropic equivalent spin chainsE3 to the case of inequivalent and anisotropic chains. 
An open question is whether in the case of inequivalent chains bound states of these elementary excitations can be 
obtained in contrast to the case of equivalent chains. A necessary condition is that there exists attractive interactions 
between a and s fermions. The study of such bound states could be of interest in relation with light scattering 



experiments on NaV20s. Besides semi-classical approximation, a nothe r approximate treatment is poss: 
case where u\ — u 2 , Si = <5a-^F and K1+K2 = 1, the Hamiltonian (3.24) is the double sine Gordon mode" 54 



n In the 

4 or the 

Bukhvostov-Lipatov modelEj. The model is-known to be Bethe Ansatz solvable^ but its elementary excitations and 
S matrix have only been obtained recentlyE3 and shown to be identical to those of the double sine-Gordon model on 
its integrable line. The spectrum . of . th &.modelEll is consists of four massive particles carrying two quantum numbers 
Qi,Q 2 = ±1, and the S matrix isESN'O: 

S = sgd®s££, (4.19) 

where S$ G is the S matrix of a sine Gordon model having the parameter K and K a = AK a /(l + 2K a ). As a result, 
one of the sine Gordon models is in the attractive regime and has a spectrum made of fermions and bound states 
of fermions whereas the second sine-Gordon model is is the repulsive regime and has a spectrum containing only 
fermions. The mass of the fundamental fermion is m = -, £-> — e^t. For K\ = 1/2, one recovers the equivalent chain 

47ra sin(7rKi ) > ^ 

and the mass of the bound states is then: 

777„ = -r— sin(7rnl£i), 77 e A, nK x < 1/2 (4.20) 
Ana 

Since for 1/4 < K\ < 3/4, the operator sin2</>i sin 202 is the most relevant, it is reasonable to expect that neither 
the marginal perturbations due to ui ^ ui nor the less relevant perturbations #1,2 ^ change the gapped nature of 
the spectrum. In this regime, one should not observe any bound state, since the 77 = 2 bound state only exists for 
K < 1/4. In terms of the original spin chain, the double sine-Gordon regime should be accessible if one chooses to have 
one chain with J z < J and the other chain with J z > J in such way that K\ + K2 — 1 and J± <C J, J z - However, the 
double sine-Gordon regime would not be observed in the spin tube case, in which we should have K\ + K 2 = 3/2. We 
expect that in the spin tube case, the system has more quantum fluctuations than in the double sine-Gordon regime. 
Therefore, no bound states of spinons will form. This , heuristic argument agrees with the numerical calculations on 



the spin tube that show no bound state of spinonsli 5 . 



V. EFFECT OF A MAGNETIC FIELD ON THE DIMERIZED PHASE 



In this section, we discuss the effect of the application of a magnetic field on the dimerized phase. In general, the 
application of a magnetic field to a gapped one dimensional spin liquid system results in the closure of the gap for a 
magnetic field of the order of magnitude of the gap. Below the critical magnetic field the magnetization is zero. Above 
the critical magnetic field the magnetization increases ras (h — he) 1 / 2 , the system is in a single component Luttinger 
liquid state and has incommensurate spin correlationsEa. Here, we have to consider two a priori different cases. In 
the first one and most academic, the two coupled chains carry real spins that couple to the magnetic field. In the 
second case, only one chain carry a spin that couple to the magnetic field, the other one carrying only a pseudospiu. 
degree of freedom that does-pot couple to a magnetic field. This last case corresponds to the spin-orbital modelsE£l 
and to the spin tube modelEl The behavior of the spin-orbital model in a magnetic field has been discussed in Refs. 
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P7| and |5^. Both references study models with SU (2) x SU (2) symmetry in the absence of the magnetic field. It is 
shown that in the vicinity of the SU(4) point, the spin-orbital model becomes equivalent to the 0(4) Gross Neveu 
model describing the gapped modes plus a c = 1 CFT that describes the gapless magnetic modes. Also, in Ref . |59| , 
the weak coupling case has been discussed. The case where both chains carry spin will be dealt with in Sec. V A and 



the case of a single chain carrying spin in Sec. VB 



A. Both chains carry spin 



If both chains carry spin, the coupling to the magnetic field (taken parallel to the z axis) is: 

' dxd x (4>i + 02 ) 



(5.1) 



In that case, we can use the decoupling of Sec. IV A , Eqs. and (4.6). The coupling to the magnetic field is 

-h 



V%— I dxd x 

7T 



(5.2) 



1. The case of equivalent chains 



We consider the case where the two chains are equivalent. This problem has been discussed in the context of spin 
ladderscd. The gap in the antisymmetric modes <p a is not affected by the presence of the magnetic field. On the other 
hand, for a sufficiently large magnetic field h > h Cl , the gap in the symmetric modes closes and the magnetization 

m = —^d x 4> s behaves as m ~ (h — h^) 1 / 2 for h > h ci . Moreover, it can be shown that the low energy modes of <p s 
are described by the following Hamiltonian: 

dx kjCfrn^ + ^W 1 



H s = 



2n 



(5.3) 



And the exponent K s takes the universal value 1/2 at h = h c . Correlation functions-in the incommensurate phase 
can be obtained in the isotropic case by a calculation similar to the spin-ladder caseE3. This time, ((j) 2 ) is finite so 
that e l0a /^ has exponentially decaying correlations. As a result, one has: 

(S+(x)S-(0)) ~ e -*/ c ,a;«£ 

A simplified expression for the operator S z = Sf + S% is: 



(5.4) 



S z = d x 4> s + cos(v / 80 s — 2ttuix) 

7T 



and the correlations of S z are: 

{T r S*(x,r)S z (0,0)) 



(urf 



constant 



x 2 + (ut) 2 



2K ' 



cos(27rma;) 



{x 2 + (ut) 2 ) 2 

There are subleading power law corrections at Q = 2irmx. As shown 
missed if one naively neglects the band curvature after refermionizatio 
by using the Haldane expansion of the spin operators and retaining the terms up to 4fcp- as we did here 



(5.5) 



(5.6) 



i,Furusaki and Zhang, these corrections aie 
il. These corrections can also be obtainedE2l 



2. The case of non- equivalent chains 

In the case of non-equivalent chains, the problem gets more complicated. The magnetic field still couples only to <j) s . 
However, when the magnetic field h exceeds the field h Cl needed to close the gap in the s modes, the appearance of a 
non-zero (d x 4>s) creates an effective magnetic field that couples to d x (f> a - If this magnetic field is not strong enough to 
close the gap in the 4> a modes, the system remains in a one-component Luttinger liquid for fields h not much stronger 
than the critical field h Cl . For h sufficiently large, the generated effective field can close the gap in the a modes 
leading to a two component Luttinger Liquid. Since no experimental ladder system with a biquadratic exchange 
much larger than the quadratic exchange much larger that the quadratic exchange and made of two non-equivalent 
chains is presently available, this two-step transition to a two component Luttinger liquid is unlikely to be observed 
exp er iment ally. 
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B. Only one chain carries spin 



This case includes the spin tube and the Mostovoy-Khomskii model and is therefore the most relevant physically. 
This case has been discussed for two coupled SU(2) x SU(2) symmetric chains in Ref. |55^. In the case where only one 
chain, say chain 1 to fix notations, carries spin, the interaction with the magnetic field is given by a term: 



- / dxd x (f>i, 

7T 



(5.7) 



In the case of equivalent chains, we can use the same decoupling of Sec. IV A , Eqs. ( |4.5| ) and (4.6) as in the preceding 
section. However, there is an important difference. Now, we have the coupling to the magnetic field in the form: 



ttV2 



dxd x ((f> s + (f>a) 



(5.8) 



As a result, now the magnetic field couples to bothEa <p s and cp a . Moreover, the strength of the couplings is exactly the 
same. We start with the discussion of equivalent chains. Then, we discuss non equivalent chain. We will show that 
in both case, the closure of the gap leads to a two-component Luttinger liquid, behavior in contrast with the usual 
spin-liquid systems that lead to a single component Luttinger liquid behaviorEa. We will also give the expression of 
correlation functions in the Luttinger liquid phase. 



1. Equivalent chains 



As a result of the symmetry between (p a and (p s , the gap in the symmetric and the antisymmetric mode close simul- 
taneously, leading to a two component Luttinger liquid ground state under strong enough magnetic field. Contrarily 
to the case where the magnetic field couples to both 1 and 2 spins, there is no intermediate single component Luttinger 
liquid phase. We expect that the magnetization m = —d x <pi/n behaves as (h ~ h^) 1 / 2 close to the threshold. It is 
also important to note that, the two sine-Gordon model being equivalent, one has for the fixed point Hamiltonian 
u* a = u* s and K* = K*. As a result, the fixed point Hamiltonian is invariant under any rotation in the (<p a , 4>s) plane. 
We can thus write the fixed point Hamiltonian as: 



H 



dx 
2^ 



u*K*{ttU) 2 



K 



■(d x <pf 



(5.9) 



Where <p — ((pi, <p2 ) and n = (Hi,^). Moreover, in the case of SU(2) x SU(2) symmetry, we have K% = K2 = 1/2 . 
It can the be shown easily, using the refermionization procedure that K* = 1/2 for any magnetic field. Equation (5.9) 
has important consequences for the correlation functions, which are of the form (for equivalent chains): 



(S+(x,t)S~ (0,0)) = (-l) x (x 2 -t 2 ) - ^ + constant 



-2l7T7n a X 



J2iirm a x 



(x - 1) 



2(^+4^7) 



((S z a (x,t) - m a )(S z a (0,0) - m a )) = cos(7ra;(l - 2m a ))(x z - t 2 )~ Ka + constant 



( x + t ) 2(K - + ^ 
1 / x 2 +t 2 



Att 2 V (x 2 - t 2 f 



where the index a indicates the spin in chain 1 or 2, m a is the magnetization, with mi = m (total magnetization) 
and m2 = . From the above expressions we deduce the following: the correlation function parallel to the field, 
(SfSf), has a staggered part shifted from the wave vector (7 = 7rto<7 = 7r(l — 2m), while the correlation function 
perpendicular to the field, < S^Si >, has an unshifted staggered mode and the uniform magnetization mode shifted 
to q = 2nm. The correlation functions for the spin of type 2, instead, are completely unaffected by the presence of 
an external magnetic field. 



2. General case 



In this section, we consider the case where the two chains are not necessarily equivalent. In particular, this is the 
case that is realized in the spin-tube under a magnetic field. 
In this case, is convenient to use the rotation, 
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n„ 



n-r<P2)/^/2, 4>a — (0i - 4>i)/V2 

(Hi + n 2 )/V2, n a = (ni - n 2 )/V2 



(5.10) 



to bring the Hamiltonian to the form: 
~ J 2n [ 2 

U\K\ - U 2 K 2 2tt tt 



(tHI,) 2 + (7rn a ) 2 ] 
ui/K-i - u 2 /K 2 



(u 1 /K 1 + u 2 /K 2 ) 



(^</> s ) 2 + {d x 4>a) 



d x <t>sd, 



xWs ^x tya 



^9 f 

dx[cosV8<f> a 



cos 



h 

71 



dxd x ((/) s + (i> a ) 



(5.11) 



Shifting (p a — > a + 7T/V8, renders Hamiltonian invariant under the interchange of (/> s and a . As a consequence, 
the gaps in <p s and a are identical and have to close simultaneously. Therefore, as in the case of identical chain, 
there is a transition from a gapped phase into a gapless two-component Luttinger Liquid phase. The most general 
Hamiltonian for a two-component Luttinger liquid is: 



H = / — [A/ aa (rfI a ) 2 + 2M as rfI a7 rrL + M ss (irll s ) 2 + 
J 27r 

+N aa {d x <t> a ) 2 + 2N as d x (f> a d x 4> s + N ss (d x s ) 2 ] . 



(5.12) 



In the case of coupled XXZ chain, the effective Hamiltonian (5.12) can be simplified by making use of symmetry. We 
know that the original Hamiltonian, Eq. ( |5.11 ) is invariant under the transformation: 



n„ 



n.. 



(5.13) 



This symmetry must be preserved by the renormalizcd Hamiltonian. Therefore, one must have M aa = M ss and 
N aa — N ss . As a result, the effective Hamiltonian (5.12) is diagonalized by returning to the original variables (/)%, <f) 2 . 
The effective Hamiltonian is therefore the following: 



H 



dx 
2^ 



Kl 



k(d x 0i) 2 + u* 2 K;(irIL2f 



K. 



Hdxhf 



(5.14) 



The absence of coupling between <p\ and 2 in this Hamiltonian is somewhat surprising. This can however be 
understood by the fact that the Hamiltonian even in the presence of a magnetic field is invariant by a rotation by 
7r S 2 ' z — * —S 2 ' z . Thus, (5f) = for any h, implying that d x (j)id^a terms cannot appear in the two component 
Lutti nger liquid Hamiltonian. In contrast to the spin ladder caseEil, the exponents K* 2 in the Hamiltonian Eq. 
(5.14) are non-universal. This can be shown in the following way. Let us consider the case where u\K\ — U2K2 and 
u\/K\ — U2/K2 are small compared to u\ + it 2 , namely where the couplings in the two chains are nearly iden tical. 
In that case, we can neglect in a first approximation the terms n a Hb and d x () 



in the Hamiltonian 5.11. We 



are thus left with two decoupled sine-Gordon models under a magnetic field . These sin&iGordon models undergo a 
commensurate-incommensurate transition at a critical magnetic field. It is well knownEil that the exponent at the 
transition assumes a universal value that renders scaling dimension of the cosine term equal to one. Therefore, in 
our case, the universal value of the exponent at the transition is: K* = 1/2. As a result, close to the transition the 
effective Hamiltonian is: 



H = 



dx 
2^ 



E 



u*K*(ttU,) 2 + — (dx^f 



dx 
27 



u l K l 



u 2 K 2 2 

7T ILjlla 



ui/Kx - U2/K2 



d x (j)sd x (j)a 



(5.15) 



Returning to <px and <f) 2 , one obtains an Hamiltonian of the form Eq. ( 5.14 ) with 



K{ = 
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indicating that except in the case of equivalent chains, xxne should not expect universal exponents at the transition. 
This should be the case in particular for the spin-tubeEj. It may provide an experimental test for the spin-o rbita l 
mode£a'c3 of NaV20s since the exponent K\ controls the temperature dependence of the NMR relaxation rateESEil. 
However, since the transition temperature between the gapped and the gapless phase in NaV2 0s is T c = 35K, the 
magnetic field needed to close the gap should be of order 52. 5T which could make the experiment impossible. In 
Na2Ti2Sb20, with T c — IIOK, the situation is even worse. 
The spin and pseudospin operators are, 

S z [x) = Tn — -^-i + — sin(2^i - 2nmx), S + {x) = ^= [e l ™ a + sin(20! - 2wmx)] 



it it a -J it a 



r z ( ;c )=-Mi + ^l sin 2^ T + (x) = -^=[e wraB +81112^2], (5.17) 

ir na - ■ 



where to is the total magnetization. n In the case of the spin-orbital model, the spin-spin correlation functions are 
therefore giv en b y the usual formulasu. The situation is however more interesting in the case of the spin tube. Using 
the formula (fO|), one has: 



(TSP(x,t)SP(0,0)) = (TS p (x,t)S fi {0,0)) x 
where (3 = (+, — , z), and p is the chain index. Explicitly, 



(Ts + (x,t)s~ (o,o)> = {-\) x / a {x 2 - (mtfy^r + const {x 2 - ( Ul t) 2 y { ^ 



i + |(Tr+(x,t)r-(0,0)) 



(5.18) 



+ 



(x — Ult) 2 (x + Ult) 2 

({S z {x, t) - m)(S z (0, 0) - to)) = cos(7rx(l - 2m))(x 2 - ( Ul t) 2 )- Kl + const ' ' 



„2 / , r^-Ki . K l ( 1 , 1 



An 2 \{x — uit) 2 (x + uit) 

(Tr+M)r- (o,o)) = (-if (x 2 - (u 2 t) 2 y*k + const { x 2 -t 2 y^ +K ^ ( (X + (Ilff ) ■ 

This correlation function, which enters in particular in the calculation of the NMR relaxation rate, contains power- 
law divergences at wave vector q ~ 0, 7r(l ± 2m) but also ±2m due to the fluctuations of chiralities. 



VI. CONCLUSIONS 



We have presented a field-theoretical analysis of the low-energy physics of the anisotropic spin-orbital model and 
the three-leg ladder with periodic boundary conditions (the spin tube) in the strong interchain coupling limit. We 
gave a derivation of the field theoretical model from the lattice Hamiltonian and then analyzed the phase-diagram 
using renormalization group equations. The system is found to exhibit a gapless phase, a spin-liquid phase or an Ising 
Antiferromagnetic phase depending on the microscopic couplings. The spin liquid ground state is two-fold degenerate, 
formed either by singlets of spins on even bonds and singlets of orbital pseudospins (chirality in the spin-tube case) on 
odd bonds or the other way round. The antiferromagnetic phase competes with the spin-liquid and we have discussed 
this competition briefly. The spin liquid phase obtains in the case of the spin tube with SU (2) symmetry, and we 
discussed the nature of excitations above this spin liquid ground state. These excitations have spin S z — ±12 and 
pseudospin r z = ±12. They are formed by introducing a free spin as a defect in the spin singlet pattern as well as a 
pseudospin in the pseudospin singlet pattern. These excitations lead to a kind of non Haldane spin liquid analogous 
to the one discussed by Nersesyan and TsvelikE£l. An interesting consequence is the absence of a magnon peak at 
q = — in the spin-spin correlation functions of the spin tube. This behavior could be tested in numerical simulations. 
We have investigated the effect of an applied magnetic field h on the dimerizcd phase. A strong enough magnetic 
field h > h Cl causes the closure of the gap and the disappearance of dimer order. The resulting gapless phase is a 
two component Luttinger liquid in contrast to the one component Luttinger Liquid that is observed in spin ladders. 
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The exponents appeared to be non-universal at the transition point, in contrast with the spin ladder casdH3. It would 
be interesting to obtain numerically the Luttinger liquid exponents for the spin tube or the anisotropic spin orbital 
model. For the spin tube, we have also shown that new soft modes appeared in the spin-spin correlation functions 
above h Cl by comparison with the soft modes of the single chain. This is the result of the presence of soft chirality 
modes. This may be tested numerically. 
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APPENDIX A: MEAN FIELD TREATMENT OF THE TWO XY CHAINS COUPLED WITH A 
BIQUADRATIC EXCHANGE RESPECTING THE XY SYMMETRY 

1. Equivalent chains 

We start from the Hamiltonian: 

H = Sn.gS n +l,a + ^n,a^n+l.a + ^ ^ (^n,2^n+l,2 + ^n^n+l^) (^n,l^n+l,l + ^n.l 1 ) (Al) 



After the usual Jordan Wigner Fermionization, this gives: 

H = a n,a a n+l,a + a n,a a n+l,a + ^ ( a n,2 a n+l,2 + a n,2 a n+l,2) ( a n, 1 a n+l,l + a n, 1 a ri+l, 1 ) 



(A2) 



The mean field approximation is obtained by taking: (recall, o = 1, 2 is a chain index). 
The mean field equation are then: 



1 = -Jl. , 

7T V V 1 + A< 



2K 1 



1 + At 



l+At 



i = — (l + AiVl+ . , , 

7T ^ V V 1 + At 



1 - 



l+At 



(A3) 



(A4) 



Where K, E a re the complete elliptic integrals of the first and second kind respectively. For small 6, the mean field 
equations (A4) reduce to : 



t 



l + At 7T 
1 



4A , / V8 I 1 + Xt 

v ln m 



(A5) 



And one obtains: 



t = -- 



7T + 2A 



S = V8- — ^ cxp f-7v 
7T + 2A V 4A 



(A6) 



This result from the mean field theory agrees with the prediction of bosonization since bosonization would predict 
the same essential singularity in 5 at small A, the interchain coupling being marginal. 



2. Non equivalent XY chains 

In this section, we consider a problem with a Hamiltonian of the form: 

H = / ](Sn^n+l + S n S' I | +1 )(1 + a(r^T n+1 + T n t^ +1 )) 
n 

After a Jordan Wigner transformation, the Hamiltonian becomes: 



(A7) 
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H = ^Z( a n a n+i + a n a+ +1 )(l + a(b+b n+1 + b n b+ +1 )) 



(A8) 



The mean field approximation is now: 

( a n a n+l + a n a n+l) = f l + H™^ 

(b+b- +1 +b-b+ +1 )=t 2 + (-r5 2 

The mean field equations are now: 



ti 



_(l + 0(ta) W 1+ _ 



■E 1 



l+crt 2 



<5i 



2a<5 2 



'1 + 



1 + at 2 



2K 1 



l + at 2 



E 1 



-( 
-(-&_Y 



^i+at 2 y 



*2 = --(i + ot 2 )yi + 



fE 1 



«2 = - 



2aJi 



-1 



2K | 1 - ( | 



E 1 



(If) 
-(*)" 



-GO' 



for small a, the mean field equations for ti and 5i reduce to: 



(A9) 
(A10) 



(All) 
(A12) 



h = - 

4aS 2 , / V8 | 1 + c^2 
oi = In — — 



(A13) 

One sees that a finite bandwidth 2a /it is produced for the r spin waves at least for small a. The gap equations admit 
solutions at small S. One expects at weak coupling an essential singularity in 5i <~ exp(— Cte/a). Therefore, 5\/t\ 
should go to zero for a — > 0. This implies that one can write: 



2a<5i 



*2 



In V8 



2crf 



nil 



(A14) 
(A15) 



One can solve the mean field equations for 5\ and 8 2 . One obtains: 

4 



* = _(! + (D) exp(-£) 
* = (l + (l))exp(-£) 



(A16) 
(A17) 



We obtain therefore self-consistently a gap much smaller than the smallest bandwidth. As a consequence, the 
correlation length is much larger than the lattice spacing, which justifies a continuum approximation. The bosonization 
treatment is therefore valid for a — > and close to the XY limit. 
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APPENDIX B: RENORMALIZATION GROUP EQUATIONS DERIVED BY MOMENTUM SHELL 

INTEGRATION 



1. The Quantum sine Gordon Model 

In this section, we derive RGE for the quantum sine Gordon model using the method of Knops and Den Oudenll. 
The quantum sine Gordon model is defined by the following lattice Hamiltonian: 

H = E t(^) 2 + - *) 2 ] - p^a cos4 ^} ( B1 ) 

The cutoff in this model is only on space and not on time. 

The Euclidean action of the Quantum Sine Gordon model is obtained as: 

S *' " £ £ ll T 5Sf (*" 4) "" 2 " H / ** » 84 * <B2 > 

Where f dx —> a^Z,-. O ne can integrate from — oo to oo using a cutoff function ip(k) = ©(— — |fc|). T-his action 
breaks rotation invariance in the x, r space in contrast with the action of the classical sine-Gordon modclcJ. Instead 
of the sharp cutoff, one can use any cutoff function <p(ak) such that ip(0) = 1 and lim^^oo <p(x) = oo. An example is 
(p(ak) = e-°l*l. j-, 

the technique developed by Knops and Den OudencZl for the classical sine-Gordon model is straightforwardly adapted 
to the Quantum Sine-Gordon model. One obtains: 

^ = (2-4^) 5 (B3) 

l(#)-(^)T** , "- MO f w ,B4) 



dp 

J W = iWT^-»l& (B5) 
at u \iru/ j op 



Where: 



Jo n 
F 2 {p)= / 2 \os 2 #e 16X [ y ^ + ^ 



IT 



(B6) 



And: 



2 

For instance, if = e~ K , one has: 



1 f°° Hk 

V(p,B) = - —(p(K)(cos(Kpcos6)e- Kp \ sine \ -I) (B7) 



v( P , o)=- 1 -J ^^uMA^m (B8 ) 



2 



The Eqs. (B3) have to be contrasted with the usual RG equations for the sine Gordon model with a cutoff 
isotropic in the x,ut space. In the latter case, the equations are such that du/dl — 0, and only K is flowing under 
renormalization. In the case we have considered, both u and K are flowing under RG transformation. However, for 
K = 1/2, the RG equations (B3) are considerably simplified. Since i*i(0) = F 2 (0) = and Fi(oo) = ^(oo) = 1, one 
has the following simplified RG equations: 
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dl \TTuJ 



(B9) 



These RG equations are identical to the ones obtained with a cutoff isotropic in x, ut space, we conclude that in 
the vicinity of the BKT transition point, the anisotropy between space and time does not matter. 



2. Renormalization of the dimerization term 



We consider the case where 5\ = 5 2 = 0, with Euclidean action: 

' u a (d x (f) a ) 2 (d T <t) a ) 2 2g 

I UXUI 

a=l,2" 



S = / dxdr 



2irK a 2-Ku a K a (27ra) : 



■ sin 20i sin 202 



Using the method of Knops and Den Ouden, we obtain the following RG equations: 



§ = <2- Kl -K, )g 



d ( ui 
d ( 1 



dl \u\K\ 



— ( H^L \ - U2 ( g 



,m r I * / dpp^- K ^d p F 5 { P ) 

dl \K 2 / 8 V 7ru 2 / Jo 



1 



. , =- f-*-Y 

dl \u 2 K 2 J Su 2 \7ru2J 



/ d PP ^- K ^d p F 6 (p) 
Jo 



Where: 



And: 



Fs{p)= I — cos 2 OexpfAId V(p, 6,1) + 4K 2 V{p,6,—) + 2(K 1 + K 2 )\np 
F 4 (p)= ( —sm 2 eexp(4K 1 V(p,e,l)+4K 2 V(p,e,—) + 2(K 1 + K 2 )lnp] 

Jo * \ Ml J 

F 5 {p)= I — cos 2 OexpilKiVfaO,—) + ^K 2 V{p, 6,1) + 2{K 1 +K 2 )\np 

JO 1 V "2 

/•2tt in / 

F 4 (p)= / —sm 2 6exp[4K 1 V(p,6,—) + 4K 2 V(p, 6,1) + 2(K 1 +K 2 )\np 

Jo ^ V u 2 

1 f°° dK 

V(p,6,a) = - / — ^(K)(cos( K pcos0)e- KQ "l sine l - 1) 

2 Jo « 



(BIO) 



(Bll) 



(B12) 



(B13) 



It is easily seen that similarly to the case in which K\ + K 2 = 2, universal coefficients appear in the RG equations. 
These universal coefficients are respectively: 



F 3 (oo) 



F 4 (oo) 



,-2tx 
JO 



cos 2 8 



cos^ 



^] sin^ 



K 2 



2tt 



d6 



cos 2 6 + 



(-) » 



sin 2 (9 



K 2 



22 



F 5 (oo) 



2 71 



cos 2 8 



(16 



cos 2 9 



fe) si 



2 „ \*i 7T 

sin I 



F 6 (oo) = / 
Jo 



sin 2 6» 



d6> 



cos z 



sin 



One can check easily that if u\ = u 2 , the RG equations (Bll) reduce to 



d ( 1 
Jl 

d ( 1 



9 



8^ 2 u 2 

^2 



9 



dl\K 2 ) 8tt 2 u 2 



which is equation (3.24) in the spintube paper for 6i = 82 = 0. 



3. The full problem 



(B14) 
(B15) 



(B16) 
(B17) 



For the full problem, the Lagrangian is: 



L 



dxdr 



28x 



2-kKx 2iru 1 K 1 (27ra) 2 
u 2 (d x fo) 2 , {d T 4>i) 2 2S 2 



COS 401 



+ 



2nK 2 

2ff 
(2-Ka) 



2nu 2 K 2 (2na) 2 
2 2 sin 201 sin 202 



COS 402 



Applying the Knops Den Ouden Method, we obtain the following Renormalization Group equations: 



d ( u\ 
Jl 

d ( 1 



dl \u\K\ 

dl \K 2 
d ( 1 



ui I g 

8 \7TUl 
1 



dpp ^- Kl -K 2)dpF3+Ui 



Jl_ 



dpp^-^d p F 2 



\ nUj i/ Jo 

2 roo -1 / c \ 2 



Bui \1ni1) J Q u\ \1tu1) J Q 



dl \u 2 K 2 
d_ /_6j_ 

dl \1TUl 

d_ 

dl \iru 2 
± (j_ 

dl \TVUl 



1 



M-) /°° dpp 2{2 - Kl ' K2) d p F 6 + — ( — J f dpp^-^d p F x 
8u 2 \nu 2 J J u 2 \nu 2 J J 

(2-AK^-U^ 2 
nui 8 \7ru1 



(2 



4if 2 )A_l r d pp 2 ^-^d p F 7 
nu 2 8 \ttu 2 J J 



Where: 



(2 - K x - K 2 ) 



1 9 ( h 



7TUl 2 7TU\ \7TUl 



poo r />co 

/ dpp 2 ^d p F 9 + ^ dpp 2 -^d p K 

Jo XU 2 J 



d9 

2^ 



4[K 1 V(p,0,^)-K2V(p,0)]-2(K 1 -K 2 )\np 



(B18) 



(B19) 
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= ^_ e i[K 1 V(p,e)-K 2 V(p,9^)]-2(K 2 -K 1 )lnp 



F 9 (p) = / e *Ki[v(P,»)+?*P\ ( B20 ) 
Jo 

One can check that if u\ = U2 = u, these equations reduce to those derived using OPE techniques. 
Clearly, the presence of a finite velocity difference results in different coefficients in the RG equations. However, 
this should not affect the topology of the phase diagram. 
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FIG. 1. Cylindrical three-leg ladder (spin-tube). The 
choice of the topology affects the strong-coupling limit. 




FIG. 2. Phase diagram of the spin-tube model. For 
K i + K 2 > 2 the system is in a gapless phase, while it 
is gapped for K\ + K2 < 2. The dashed line corresponds 
to the isotropic line with K\ — K%. The isotropic point 
{K\ — 1,K2 = 1) corresponds to two XY chains, while 
the isotropic point (K\ — 1/2,1^2 = 1/2) corresponds to 
two-equivalent spin-chains considered by Nersesyan and Tsve- 
lik. 
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FIG. 3. Representation of the two-fold degenerate 
ground-state, formed by singlet of spins on even bonds and 
singlet of chiralities on odd bonds or singlet of spins on odd 
bonds and singlet of chiralities on even bonds. 
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FIG. 4. Phase diagram in the limit K\ — K2 = as a 
function of g for fixed 8\ < S 2 . An intermediate phase with 
dimer order in chain 1 and mixed dimer and antiferromagnetic 
order in chain 2 is obtained when 4^/8162 < g < 4^2 • 
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FIG. 5. (a) magnetic soliton. These solitons are asso- 
ciated with the triplet excitations having m = ±1. On the 
figure, m = — 1; (b) non- magnetic soliton. These solitons are 
associated with the triplet excitation having m = or the 
singlet excitation. 
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FIG. 6. (a) triplet excitation of the spin ladder with bi- 
quadratic exchange; (b) singlet excitation of the spin ladder 
with biquadratic exchange 
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FIG. 7. Physical picture of the spin excitations above the 
ground state. The elementary excitations confine to form sin- 
glets or triplets between one unpaired spin and one unpaired 
chirality. 
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